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INTRODUCTION 


Current concern with the affects of transport on atmospheric chemistry, in particular 
the ozone chemistry, have regenerated and strengthened interest in the mass exchange 
between the troposphere and stratosphere. It is well known from aircraft experiments 
(Danielsen, 1964; 1968), synoptic analyses (Danielsen, 1959; Reiter, 1963; Reiter 
and Mahlman, 1965; Mahlman, 1965 and others) and numerical model simulations 
(Mahlman and Moxim, 1978) that stratospheric air is intermittently transported into 
the mean troposphere by wave motions during cyclogenesis in extratropical latitudes. 
Although numerical estimates of the amount of mass leaving the stratosphere each 
year by this process vary considerably, there is general agreement that it is a 
significant fraction of the total stratospheric mass. There also is general agreement 
that the corresponding inflow of tropospheric air to the stratosphere occurs at 
lower latitudes. But it remains to be determined what percentage of the inflow is at 
tropical latitudes, and whether the physical process is predominately a zonal mean 
transport by the upward branch of the Hadley cells, or a combination of mean and wave 
or eddy transports. 

The NASA -Ames Research Center's Inter- Tropical Convergence Zone (ITCZ) Experi- 
ment was organized to study tropospheric-stratospheric transport associated with 
the ITCZ by means of aircraft observations of several important trace species, in- 
cluding ozone, nitric oxide, water vapor, chlorofluromethanes and aerosols. Although 
the emphasis was on determining the mean transport rates, the experiment was 
designed to resolve spatial and temporal variations in some of these trace species 
and, in the meteorological variables, to objectively determine their existence and to 
reliably establish their frequencies, amplitudes and relative phase relationships. In 
support of this secondary objective the recording rates for several U-2 measurements 
were increased and special arrangements were made to permit high resolution deter- 
minations of the vertical profiles of temperature, relative humidity, ozone and the 
horizontal wind velocity. 

This report deals specifically with the temperature, relative humidity and wind velocity 
profiles resulting from objective computer processing methods applied to digitally 
recorded radio -ozone sonde data. In particular, the wind velocities are derived 
directly from special radar tracking data, which eliminate the major sources of 
errors that limit accuracy of the conventional rawinsonde system. In addition, by 
having a completely independent set of tracking measurements the accuracy of the 
conventional measurements can be assesed. The accuracy of the computed winds can 
be directly determined, of course, but also from differences in the radiosonde computed 
heights and the radar heights, the less direct effects of systematic temperature and 
pressure errors can be estimated. These results are of major significance in the 
tropics because the height variations are small and the signal -to -noise ratio is 
easily swamped by the instrumental and/or data processing errors. 
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The experiment involved the cooperation of several governmental agencies whose tech- 
nicians were asked to provide non-conventional data and/or special recording methods. 

It follows quite logically that difficulties would arise in the data recording and 
processing. Indeed, many were encountered. However, by persistent efforts and 
subtle detective work we overcame the obstacles and are proud of the final data set. 
There are unavoidable gaps in the processed data but the original objective has been 
met. Also, since the entire set of observations were digitally recorded at 1 or 2 
per second, it was possible to design and apply low-pass digital filters to each variables’ 
vertical profile. Most importantly, by subtracting these filtered profiles, i. e. , by 
effectively applying narrow band-pass filters, individual wave modes can be isolated. 
Thus, for the first time, we can isolate examples of the Kelvin wave mode, the Rossby- 
gravity modes, etc. , and compare their vertical wavelengths, amplitudes and phases 
to those dervied from perturbation theory. 

Recognizing that meteorologists usually work only with processed data, we are pre- 
senting sufficient background and discussion to enable one to appreciate the unique 
advantages of this data set, and how these objective processing methods differ from 
conventional methods. However, to avoid a myopic preoccupation with details, many 
explicit steps have been eliminated. They will be supplied upon request. 

RATIONALE FOR RADAR TRACKING AND DIGITAL RECORDING 

In principal, radar tracking of a radiosonde balloon system has several advantages over 
the usual GMD-1 rawinsonde tracking. Because the radar transmits a short pulse of 
energy and measures the time for this pulse to go to and return from the target, it 
obtains a direct measurement of the slant range with minimum interference from re- 
fracted or reflected waves arriving later due to their increased path. The radar always 
points toward the most direct path to the target, slightly curved by a gradient in the 
index of refraction. On the other hand, the GMD antenna can err in its direction angles 
because it senses only the direction of maximum radio signal intensity, and this max- 
imum can be the result of interference between a most direct ray and a more refracted 
or reflected ray. As shown by Danielsen and Duquet (1967), the seriousness of this 
interference error increases rapidly as the elevation angles decrease to 15 degrees or 
less. During the ITCZ Experiment the elevation angles remained larger, so inter- 
ference did not produce serious errors, but the effects of interference can be demon- 
strated. 

Another advantage of radar tracking (for this experiment the major advantage) is that 
the radar alone supplies a necessary and sufficient set of measurements for wind com- 
putations. Knowing the azimuth and elevation angles, plus the slant range from the 
radar to the target, one can directly compute the three-dimensional position of the 
target over the curved earth as a function of time. Then, by finite differencing the 
horizontal positions for a constant time interval, one can compute a time-averaged 
velocity. The difficulty with the GMD-1 is its lack of a slant range measurement. 

To determine the radiosonde transmitter's three-dimensional position over the curved 
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earth, it is necessary to compute the height from the radiosonde observations of 
temperature, relative humidity and pressure, and then use this height in combination 
with the measured azimuth and elevation angles., As will be shown, errors in pressure 
and temperature (especially errors in pressure) contribute to Large height errors at 
high altitudes where the ambient pressure is low. If the measured pressure is too 
low, the apparant slant range is too large and the horizontal distance out is also too 
large. Since an error in the computed winds depends on relative errors in position, 
an increasing error in position is most serious. Such is the case when the pressure 
error is constant; then the percentage error increases as the balloon ascends to lower 
pressures. 

In addition to the errors in the computed wind vectors, an error in pressure produces 
an error in the pressure or height to which the wind is assigned. Again, such phase 
errors tend to increase with decreasing pressure, increasing height. It follows that 
the geopotential heights and winds dervied from radios onde-GMD-1 data are most 
accurate at low elevations and their accuracy degrades as the balloon ascends. For 
these reasons it has been standard practice to increase the finite differencing interval 
above the 7 km level. For example, at operational rawinsonde stations, the azimuth 
and elevation angles are recorded only once per minute. The successive positions are 
plotted graphically, and a two-minute average wind velocity is computed by central 
finite differencing. Above 7 km this difference interval is doubled, yielding a four- 
minute average wind. Because the balloon ascends at an average velocity of 5 m sec^, 
these time averaged winds correspond to a 600 m average below 7 km and 1.2 km 
average above 7 km. 

Here, in processing the winds from the radar tracking data, we have used a constant 
200 m average, for the following reasons. Based on precision radar tracking data 
from an FPS-16 radar digitally recorded at 20 times per second, Danielsen and Duquet 
(1967) showed that a 2 sec averaged wind resolved even the torroidal motions of the 
balloon. Power spectra analyses indicated these balloon- induced oscillations had 
periods ranging from 6 to 20 sec, perhaps extending to 30 or 40 sec. In general, 
there was a distinct minimum in the power or variance spectra from 40 to 60 seconds; 
therefore, to eliminate spurious oscillations caused by eddy shedding from the 
ascending balloon, it is convenient to use a 40-second difference in the wind computa- 
tions. At an average ascent rate of 5 m sec -1 this corresponds to a 200 m difference 
in height. We changed the independent variable from time to height when we discovered 
that there were rather large systemmatic differences in the ITCZ Experiment’s balloon 
ascent rates. These differences altered the effective vertical wavelengths of the 
digital filters, an alteration which is unacceptable. Ordering the data at a constant 
5 m interval, by interpolating in the 1 second data, eliminated this balloon-dependent 
effect. Thus, the filters pass only those vertical wavelengths exceeding a prescribed 
wavelength. 

Clearly, since the winds are derived from differencing position coordinates, relative 
errors in horizontal position are more important than absolute errors. However, to 
take advantage of the additional radar tracking capability, i. e. , to use the radar 
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derived heights to assess errors in the radiosonde system, we must consider the 
absolute errors. The heights for the radiosonde are computed from the hydrostatic 
approximation and the equation of state for a mixture of ideal gases. Explicitly, the 
appropriate equation is 

/ Po 
R 

g T y d< ln P) 

r 

where Z is the height at the pressure p, R is the specific gas constant for the mixture 
of diatomic molecules in dry air, g is the acceleration of gravity and T y = T(1 + . 612 Xy) 
is the virtual temperature, augmented by the specific humidity Xy to account for the 
reduction in mean density by the lighter-than-air water molecules. Here, Xy is ex- 
pressed in grams of water per kilogram of moist air. 

The appropriate value for R is 2. 87x10 6 ergs gm'-^deg -1 , a constant. It is customary, 
also, in meteorology to treat g as a constant, 980 cm sec - ^ and then to refer to Z as 
a "geopotential’’ height. On the other hand, the height determined from the radar is 
the geometric height. Thus, before direct comparisons are made the geometric must 
be converted to geopotential height, or vice-versa. Here we have made the former 
conversion and will then assess the errors in the geopotential heights in terms of 
errors in T v and P. The full potentials of this method could not be realized in the 
ITCZ Experiment due to the limitations of the tracking radars and the difficulties the 
radar technicians encountered getting the radars on target for automatic tracking,. 

During the ITCZ Experiment we were fortunate to have the cooperation and support of 
the U. S. Army Signal Corps, who provided radar tracking of the ozonesonde balloons 
and the U-2 and Lear Jet aircraft,, Two NIKE radars, normally used for rocket 
tracking, were made available for the experiment. However, they did not have a 
digital readout capability, so technicians from NASA- Wallops Island facility supplied 
portable digital recorders and made the necessary modifications and adaptations on 
site. Considering the inherent difficulties in this task and the time limitations they 
were very successful, but occasionaly the digitizers failed, producing some gaps in the 
tracking data. 

Major delays in processing the data were caused by systematic errors in the radar’s 
tracking measurements. Errors were detected in both angle measurements and in the 
slant range. Fortunately, by comparing the GMD-1 and radar tracking of elevation 
and azimuth angles for each balloon ascent over the complete set of ascents, the relative 
errors were determined as a function of time. Then from a few cases, where both 
radars were simultaneously tracking the same balloon or the same aircraft, reasonable 
absolute errors could be assigned to both radars and to the GMD. The most troublesome 
and least understood errors were in the radar’s measurements of the slant range. The 
simultaneous tracking data sets indicated unusually large errors in slant range which 
remained constant during each flight but which differed from flight to flight. To over- 
come this problem we were forced to reverse our logic and use the radiosonde-determined 
heights to determine the error in slant range for each balloon ascent. 
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The effects of this reversal in logic are not as bad as they would appear to be. We 
used only the first minute of radar tracking data, after we were sure the radar was on 
target, to calibrate the radar slant range. Because it took some 6 to 10 minutes to 
get the radar on target, the calibration minute corresponds to pressures of 700-600 mb. 

In this pressure range the radiosonde heights are quite accurate with standard errors of 
about 10 to 20 m. As will be shown, the errors increase with height to hundreds, even 
to thousands of meters, so the initial calibration does not preclude the higher altitude 
error determinations. 

Finally, to provide tracking data for the entire ascent and, therefore, to permit com- 
putations of the complete wind velocity profiles, we transformed the GMD tracking 
angles to the radar site, eliminating the parallax errors, and converted the radiosonde 
heights to an equivalent slant range from the time of balloon release until the calibra- 
tion minute. The justification for this extension of the data comes from the excellent 
correspondence between the GMD-1 and radar tracking measurements at the generally 
large elevation angles which occurred during the experiment. Before presenting and 
discussing the results we will breifly outline the logic used to prepare and process the 
data, including the relevant equations. 

DATA PREPARATION 

The radar tracking data were digitally recorded at a rate of 2 per second, but since two 
radars were usually in operation the data from the two were alternately recorded on one 
channel. This alternation reduced the data rate for most ascents to 1 per second, so 
this rate was accepted as standard. Data preparation included correcting for system - 
matic errors in azimuth and elevation angles and a special correction for slant range, 
which was determined by the computer after the radiosonde data were prepared and 
processed. In view of the dependence on the radiosonde data, we will defer the slant 
range calibration discussion and focus first on the radiosonde data preparation. 

Because of the time constraints, the most practical operational procedure was to attach 
the ozone instrument to the balloon train and to use the radiosonde's radio transmitter 
for both sets of measurements. This means that 8 different signals are transmitted 
in 4 to 16 second bursts, in a semiorganized sequence, and all data could be recorded 
on one channel. Obviously, the system was not designed for digital recording and 
computer processing. Nevertheless, we were stuck with these limitations or no digital 
recording at all. We accepted the challenge to unscramble the data. 

The first step was to develop a logic for separating the ozone data from the radiosonde 
data. In this discussion the following symbols will be used for the 8 separate frequencies. 

OZONE DATA RADIOSONDE DATA 

OZZ - Ozone zero RTF - Radiosonde temperature 

OZC - Ozone calibration constant RHF - Radiosonde humidity 

OZB - Box temperature RLR - Low reference 

OZF - Ozone signal RHR - High reference 

Under normal operations the radiosonde data are transmitted sequentially in an order 
determined by a mechanical baroswitch. There are some 150 switches from RTF to RHF 
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or from RTF to RLR or RHR in the 1^ to 2 hours it takes for the balloon to reach its 
maximum altitude, or about 36 seconds per contact. Usually the RTF is transmitted 
for six -tenths of this time and the other for four-tenths. Again, it must be stressed 
that the system is designed for manual -graphical readout and processing which accounts 
for the long transmission times for each signal. But, of course, when only one signal 
is transmitted at a time there are large data gaps in each signal's trace. In graphical 
processing, the operator visually ignores these gaps by following the trends of the 
recorded segments. Here, we must fill the gaps after the data are separated. 

Adding the transmission requirement of the ozone instruments increases the number 
and, to a lesser degree, the duration of data gaps. Transmission of the ozone data is 
controlled by a battery -driven commutator which interrupts and overrides the radio- 
sonde's baroswitching. A complete cycle takes about 128 seconds and includes 9 four- 
second bursts of ozone data: 1-OZZ, 1-OZC, 1-OZB and 6-OZF's with the intervening 
radiosonde data bursts ranging from 8 to 16 seconds. Being battery-driven and mechanical, 
the cycle times were not constant; therefore, the computer algorithm had to include a 
method for determining the individual cycle periods. 

The method for determining periods is based on filtering the data with three compara- 
tively narrow frequency gates, which pass only the slowly -varying OZZ, OZC and OZB 
signals, followed by a systemmatic contraction and expansion of the time gates to 
minimize the number of nonzero entries from those three ozone signals. All frequencies 
recorded within the OZF’s time gates are then stored without applying an OZF frequency 
gate since this signal can vary more rap idly „ The remaining frequencies are then con- 
sidered to be radiosonde data and the problem remains to separate the four types of 
signals. 

First, narrow frequency gates are derived to isolate the RLR and RHR, reference 
signals. By referring to the ozone time gates, all gaps in RLR and RHR are filled by 
linear interpolation. Then using the prescribed sequence of alteration between RLR and 
RHR the individual contacts are numbered and checked against initial and final contact 
information coded from visual inspection of the recorder trace. From contact 5 to 
contact 135 in multiples of 5, these reference frequencies are preceded and followed by 
RTF unless the ozone data interfere. Using the means of these RTF signals, the trends 
in the temperature frequencies are determined and these trends are used to establish 
frequency gates to isolate RTF. Then all remaining frequencies are stored in RHF. 

For all contacts greater than 135, humidity frequencies are no longer transmitted and 
the transmission alternates between temperature and reference frequencies which are 
easily separated. 

Next, the RHR's, which are transmitted only to aid the operator visually identify certain 
contact numbers, are converted to RLR’s, and then the contacts which are not divisible 
by 5, between the first and No. 135, are identified and numbered. The identification is 
based on a change from RTF to RHF 0 Where an ozone signal interiors with this change 
the contact is assigned to the center of the four-second interruption. 
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Finally, the reason for identifying and numbering the contacts is to determine the pressure 
versus time. From a pressure calibration chart, the pressure is known only when each 
contact is actuated. No other pressure data are available. To complete the data prepara- 
tion all of the gaps must be filled and spurious signals eliminated. Spurious signals 
can occur if the counter is active while the frequency is switching, and when the signals 
are weak or noisy. 

The methods used to identify spurious signals and to fill data gaps depends upon the ex- 
pected variability of the signal and the number of data points in a given data burst or 
interval,. For slowly -varying frequencies, such as OZZ, OZB, OZC and RLR, the mean 
and standard deviation a are computed for each interval. Points whose deviations exceed 
3a are rejected, new means are computed and linear interpolation is used to fill the gaps 
between the means,, For more rapidly-varying frequencies, such as OZF, RTF and RHF, 
a linear root-mean-square, RMS, fit and a are computed if the data contain less than or 
equal to 8 data points (4 seconds). Again, points whose deviations exceed So are re- 
jected, the linear RMS fit is recomputed and then the value of the central frequency and 
its derivative are used to interpolate across the gap with a cubic interpolation function. 
Thus, f and 6f /6t are continuous from data burst to data burst. 

The remaining gaps in RTF and RHF, which now exceed 8 data points, are first fitted 
by a cubic RMS function. Next, the value of f and fcf /at at 0. 5 sec from the data 
intervals' boundaries are used to provide continuity for a cubic interpolation formula. 
Finally, the pressures are interpolated from the contact versus time data by linear ' 
interpolation in lnp„ This completes the ozone -radiosonde data preparation which yields 
a 10 -column matrix including: TIME* INDEX* OZZ* OZF* OZB° OZC* RTF* RHF* PRESSURE* 

The methodology for this data preparation was written by E„ Danielsen and programmed 
by Steve Gaines. Steve is commended for his careful and painstaking efforts which 
successfully unscrambled what first appeared to be a hopeless maze of data. 

Two consecutive examples of the computer -gene rated plots of the separated and inter- 
polated data are presented in Figures 1 and 2. Note that decreasing RTF corresponds 
to decreasing temperature, but decreasing RHF to increasing relative humidities. 

In these figures the slow but persistent drifts in the calibration frequencies are 
evident as is the large variability in RHF. Note the general increase in relative 
humidity during the six hours implied by the decrease in RHF's. These two figures 
illustrate also the large differences in ascent rates (compare the elapsed time to 
contact 150 or 155). 

At this point we can return to the radar data preparation,. The difficulties here are to 
identify and reject spurious data, to interpolate across data gaps when the digitizer 
failed and to determine the systemmatic tracking errors. That the latter existed and 
were large was evident when the radar and GMD tracking angles were plotted as a 
function of time on the same diagram. In Figure 3, for example, the two traces are 
almost parallel while the elevation angle varied from approximately 30 to 70°, but 
are systemmatically offset. The difference angle, plotted on an expanded scale to the 
right, includes a systemmatic offset of approximately 0. 7° which increases to approxi- 
mately I. 0° due to parallax. Superimposed cn these combined effects are high frequency 
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Figures 1 (left) and 2 (right). Ozonesonde and Radiosonde frequencies separated 
and interpolated. See text for symbols. 


oscillations which come from the GMD tracking data. These are the fluctuations pro- 
duced by interference between direct and refracted (or reflected) rays, referred to 
earlier. Notice that their amplitudes are less than or equal to 0. 2° in this example. 

Systemmatic errors in azimuth angles can be identified in Figure 4. Here the GMD trace 
in the left diagram has been offset by 2° to provide a visual comparison of relative 
tracking accuracy,, Again, the relative tracking capability is excellent. All significant 
variations are reproduced. The actual angle differences plotted to the right include a 
systemmatic difference of approximately 0„ 2° which varies from +0. 5 to -0. 5° due to 
parallax. From these data one could deduce that the GMD was located to the NNE of the 
radar. Comparing the right diagrams of Figures 3 and 4, it is clear that the interference 
effects are not as large nor as organized in the azimuth as they are in the elevation 
angles. 

To determine whether the systemmatic differences in azimuth and deviation angles were 
constant through the experiment, separate plots like Figures 3 and 4 were made for 
each balloon flight. From these plots it was determined that the calibration of the 
northern radar, hereafter referred to as Radar No„ 1, was changed twice. After the 
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Figures 3 (left) and 4 (right). Tracking angle comparisons between GMD-1 and 
Radar 2. 


first flight on the 16th the azimuth angle was adjusted by 3. 3°, and after the second flight 
on the 18th the elevation angle was adjusted by 0. 66° to agree with Radar No. 2. No 
adjustements were made to Radar No. 2 and this radar was finally judged to be most 
accurately calibrated. The GMD calibration for elevation angles was changed on the 
30th after it was recognized that a sudden jump of approximately 1° had occurred during 
the first flight. The effect of this jump is evident in Figure 5 where the elevation angle 
difference between the GMD and Radar 1 is approximately 1. 6° 0 For this flight only, 
both radars and the GMD were tracking the radiosonde target. In the left diagram. 

Radar 2’s data were offset by -1° to permit visual comparison. Without this displace- 
ment Radar l's trace overlaps that of Radar 2. Consequently, the difference plot GMD- 
Radar 2 would also overlap the righthand plot, These simultaneous measurements 
confirm the conclusion that the oscillations are in the GMD trace, not in the radar 
tracking. Again, we see the amplitudes are less than or equal to 0.2°. 

Having established the relative errors between the three systems, it was still necessary 
to assess the absolute errors. Obviously only one set of errors had to be assessed 
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TRACKING COMPARISON 30/7/79 3R0 ASCENT 


and then all others could be computed 
from the relative errors. The final 
judgement was based on the stability 
of the Radar 2 measurements and the 
fact that the GMD could be calibrated 
only to an accuracy of ^2°. The 
following corrections were therefore 
added to all measurements prior to 
data processing: 

Elevation Angle Corrections 

Radar 2 A« = -0. 35 
Radar 1 A 6 = +1. 03 thru 18th 
Ae = -0.375 19th-31st 
GMD At = -1 0 15° thru 29th 

A e = -2. 125 on 2nd & 3rd 
flights on 30th 
A « .= -0. 14 31 July 

Azimuth Angle Corrections 


Radar 2 
Radar 2 


GMD 


A a = 0 

Aa = -2 0 5°for 1st & 2nd 
flight on 16th 

A a = +0. 823 for all remaining 
flights 
Aa = +o. 20 



GMD -RRDBR 1 


Figure 5. Tracking angle comparisons 
between GMD-1 and Radars 1 and 2. 


Finally it was necessary to determine the accuracy of the slant range measurements, and 
this was not an easy task. Initially we were told that there were no slant range errors. 
However, when the heights computed from the radar were compared to those computed 
from the radiosonde, it was obvious that the differences greatly exceeded the errors in 
the radiosonde. Communications with the radar technician produced systemmatic errors 
which had inadvertently been overlooked. Further work revealed that these errors were 
off by an order of magnitude, but still the data were inconsistent from flight to flight. 

At this point we turned to the aircraft tracking data, selecting all examples when both 
radars tracked the same aircraft. On the basis of these simultaneous tracking sets it 
was determined that the slant range measurements were inconsistent but that the 
relative range error remained almost constant during each flight. Most of these tests 
involved low elevation angles, therefore, the distances -out were not sensitive to errors 
in the elevation angles. 

Because the slant range errors held constant during each flight but varied from flight to 
flight, we decided to use the first minute of reliable radar tracking data to calibrate the 
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slant range for each balloon ascent. Reference to Figures 3 to 5 indicates that approxi- 
mately ten minutes elapse before the radar is on automatic tracking. For the earlier 
ascents, approximately 20 minutes were required by the radar operators to find the 
target - so the calibration minute is always at altitudes where the radiosonde accuracy 
is of the order of ^Om. All subsequent height errors exceeding these limits can be 
determined. 

DATA PROCESSING: RADIOSONDE AND RADAR DATA 


The radiosonde temperature and relative humidity frequencies are first corrected for 
frequency drift, determined by the change with time of the low reference frequency. 

Then the temperatures and relative humidities are computed from the equations 
supplied by the radiosonde manufacturer. Details of the method will be supplied upon 
request. Having determined the temperature, the computer uses the Smithsonian 
Tables water vapor equations to determine the saturation pressure over a water 
surface. Multiplication by the relative humidity yields the actual vapor pressure. 

Then the specific humidity Xy is computed and used to calculate the virtual temperature. 

T v = T(1 + .612X y ) (2) 

where T is expressed in °Kelvin andXy in grams kgrams - -*. 

The following potential temperatures are computed also: 

1. Potential Virtual Temperature 0 y 



where p is the total pressure. 


2. Equivalent Potential Virtual Temperature 9 ve 


9 


ve 



(4) 


where T c is the condensation temperature, the latent heat L is a function of T c and c 
is the specific heat at constant pressure. A special numerical method for computing** 
T c was developed for these computations which yields excellent accuracy in two steps. 
Convective instability, so important in the tropics, is determined from the vertical 
derivative of 9 ye . Instability is present when this derivative is negative. 


The geopotential heights, of course, are computed from (1). Because the data are at a 
constant 0. 5 second interval, a linear mean of T y is sufficiently accurate for each 
integration step. Again we stress that these heights are geopotential not geometric 
heights. 

The radar data require special treatment. Occasionally, when the digitizer failed, 
zeros were recorded for azimuth, elevation and slant range, or values somewhere between 
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zero and the proper signals. To identify these, the second differences are first computed 
and then the standard deviation of these second differences. Values exceeding 3o are 
then examined for successive sign changes. When data gaps are identified they are 
filled by linear interpolation. In a few cases exceptionally noisy data segments proved 
to be too difficult to correct by objective computer methods, so special manual adjust- 
ments were necessary. All corrections are indicated in the printout. 

After the radar azimuth and elevation corrections are made the slant range correction 
is evaluated. For each second of the 1 minute calibration segment the geopotential 
heights of the radiosonde are converted to geometric heights and these heights and 
elevation angles are used to compute the slant range from 

r = -RgSinC +\/(R e sin6) 2 + (Z 2 + 2R e Z) (5) 

where R e is the radius of the earth and Z 
is the geometric height above sea level 
(see Figure 6). The observed radar slant 
ranges are subtracted from these com- 
puted values and the mean value is used 
to correct all radar slant range measure- 
ments for that particular balloon ascent. 

From the time of the balloon release 
until the time of the calibration minute 
the missing radar measurements are 
supplied from the radiosonde measure- 
ments. To assure continuity in the 
measurements at the interfacial time, 
the GMD tracking angles are calibrated 
by transforming the radar angles, during 
the minute overlap, to the GMD site 
removing the parallax effects. Sub- 
tracting these transformed angles from 
the GMD angles provides a correction 
for the latter. Generally, these 
corrections are small, less than or 
equal to 0. 2°. Again, the values used 

are printed out for each ascent. After these corrections are determined the radiosonde 
heights and elevation angles are used to compute slant ranges, then these and the 
azimuth and elevation angles are transformed to the radar site to complete the radar 
data set. 

At this stage the data are ready for filtering with time as the independent variable. The 
filters to be discussed in the next section were initially designed as low -pass filters 
with their 50% transmission at 0. 5, 3 and 9 minutes. For the reason discussed earlier 
we changed from time to height as the independent variable. Thus, instead of ordering 



Figure 6. Slant range r as a function of 
height Z and elevation angle e . R 0 is 
the radius of the earth. 
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the data, including the height at 1 second intervals, the data and time are ordered at 
5 m intervals,. Having established a complete set it is a simple matter for the computer 
to determine the geometric height from the elevation angle and slant range 


Z — R. 


i 1 + R 


2r sin6 + r 2 


2 -1 


and the distance out over the curved earth 


( 6 ) 


D = R e 9 = R 


cos 


_1 /cos 6 


1 +- 




-e 


where 0 is defined in Figure 6. 
FILTER DESIGN AND APPLICATIONS 


In keeping with our objective to detect wave modes in the ITCZ data, and to be able to 
study their amplitudes and phase relations in the height rather than the vertical wave 
number domain, we have designed several low-pass, non- recursive filters and 
applied them to each of the computed variables. To significantly reduce aQy aliasing, 
the response of each filter is designed to closely approximate the desired response 
illustrated in Figure 7„ Here k = 2x/L where L is the vertical wavelength and Lq is 
the grid distance or vertical distance between successive data points. For our filters 
Ld is greater than or equal to 5 m, the min imum grid distance. 



K,-K — K 0 


Figure 7„ Response function for digital low-pass filter. 


Each filter will be identified by the wavelength associated with the 50% response which 
corresponds to = -Q + = 2?r 

2 l 1/2 
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Applying two filters with k]y 2 (l) greater than ^i/ 2 ( 2 ) and subtracting the output of the 
second from the first is equivalent to applying a band-pass filter as illustrated in Figure 8. 



Figure 8. Band-pass filter derived from two low-pass filters. 


The width of this band-pass depends, of course, on the difference between k-j/ 2 (l) and 
ki/ 2 < 2 ) and upon the width of the transition zone kg-k^ for each filter. As will be shown, 
reducing the width of the transition zones tends to increase aliasing because the actual 
response for a filter with finite number of weights oscillates about the desired response. 
To reduce the aliasing it is then necessary to increase the number of weights, which 
increases the distance required for the filter. 


When the filter is applied to 2N+1 data points the filtered value at height index Jq is given 
by 

N 

X(j 0 ) =f(0)XG 0 ) + D f(j')[xo 0 +j') +XG 0 -j')l (8) 

r=i 

N 

where j r =j-jQ. The weights f(j), normalized by ^ f(j) = 1, are derived from the Fourier 
transform of g(k), see Figure 7. j=-N 
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The first bracketed term is the Fourier transform of a step function, the step being 
located at k = (k 0 -Hc^)/2 = k 1 / 2 » see the dotted line in Figure 7. The transform 
(sin k-jy 2 Z)/Z represents a damped oscillation with its maximum amplitude at Z = 0. 
The weights corresponding to this term will oscillate in sign and damp in magnitude 
away from the central point of the filter. As ^i/2 increases, the relative importance 
of the central weight and its closest neighbors also increases, so the total number of 
weights can be decreased accordingly. However, when k-^ is small, as in most of 
our applications, excessively large numbers of weights are required for such a step- 
type filter. 

The second bracketed term determines the effects of the finite transition zone. As 
kQ-k-^ increases, the square term in the denominator eventually dominates and then 
rapidly reduces the magnitude of all subsequent weights. Therefore, increasing the 
width of the transition zone permits a reduction in the number of weights for a given 
threshold of aliasing, but it decreases the slope of the filter response at k^/2* 
Generally, one must compromise and increase kg-k^ as k ^/2 decreases to avoid use 
of excessive numbers of weights. On the other hand, if the number of weights is 
kept constant as kjy 2 decreases then kj-j must also decrease, i„ e. , the data interval 
increases by reducing the density of points,, The effect is much the same as before 
but fewer computations are required. 


Changing the independent variable from height Z to the index j eliminates the height 
interval Lq in all terms except the linear term in the denominator of the first bracket, 
thus (9) becomes 



(10) 


where M is the number of data (height) intervals between L-jy^and Lq, and T is the 
number of intervals between the effective transition wavelength Lrp and Lq. This 
effective wavelength is defined by Lrp = LqL^/L-^-Lq and it approaches infinity when L^ 

and Lq ►- L -|y 2 . From (10) we can see that the normalized f(j) will be independent 

of Lq. Therefore, the weights for a given filter depend on M, T and N. If these three 
quantities are held constant, increasing l d increases the wavelength at the 50% response 
and spreads the transition zone over a larger wavelength range. 

For the first three filters illustrated in Figure 9 the data interval L = 5 m is kept constant 
and the number of weights n = 2N+1 is increased a k^y 2 is decreased. Filter No. 1, 
which passes the largest wave numbers, shortest vertical wavelengths, has M = 30, 

T = 37. 5 and N = 90. Filter No. 2 has M = 180, T = 140. 6 and N = 162. Filter No. 3 
has M = 540, T = 500 and N = 453. Since Lq = 5 m, the wavelengths corresponding to 
the 50% response are L-jy 2 (l) - 150 m, L-^^ 2 ) = 900 m an d l i/ 2 ( 3 ) = 2. 7 km. At an 
average ascent rate of 5 m sec - * these correspond to 30 second, 3 minute and 9 
minute filters. These filters are applied to the primary data sets, then the secondary 
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Figure 9. Actual filter response for the low-pass filters used in 

data processing. 

data are computed from the filtered values. For example, the winds are computed from 
the finite differences between filtered positions and filtered times corresponding to *100 m. 
This means that all reported winds are finite differences over 200 m or approximately 
40 seconds. 

After studying many of the plotted profiles of the filtered wind components, including the 
band-pass profiles, it seemed desirable to extend the filtering to still longer vertical 
wavelengths. In particular, the midpoint wavelengths were increased to 8. 1 and 24. 3 km, 
to determine whether the Kelvin wave could be separated, and to determine its vertical 
wavelength. Decreasing k-jy 2 to such low values increased M from 30 to 1620 and 4860 
for Filters 4 and 5. To avoid using correspondingly large numbers of weights we reduced 
M by a factor of 20, by increasing L D to 100 m, and then limited N to less than or 
equal to 149. Despite the rather poor filter response of the 8. 1 and 24. 3 km filters 
(see Figure 9), the results were rewarding and we will discuss them briefly,, Also, 
they stimulated further refinements of the long wavelength filters, including the intro- 
duction of another filter with the 50% response at 13 km. The response for these three 
filters, Nos. 4’, 5' and 6', are plotted as dashed lines in Figure 9. Being much steeper 
than Filters 4 and 5, they also oscillate more around 0 and 1. 
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Before discussing the results we want to credit Steve Hipskind for an excellent job of 
writing the computer codes for processing the radiosonde and radar data and for writing 
all of the plot codes. 

PRESENTATION AND DISCUSSION OF RESULTS 


A total of 40 soundings were processed 
and plotted for this report. As indicated 
earlier, failures of the digitizers were 
responsible for most of the omissions. 

A few cases were omitted because of 
erratic radar slant range measurements 
which produced erroneously large winds 
and balloon ascent rates. Using the 
third sounding on 24 July as an example. 

Table 1 reproduces the first 50 lines 
of the computer output from Filters 
1 to 3. The results are printed at a 
constant height interval of 100 m. These 
results are stored on a disk at Ames 
Research Center and will be made 
available upon request. Note that the 
last column gives the ascent rate of 
the balloon, averaged as are the winds 
over a 200 m interval. The fourth 
column gives the geopotential height 
error between the radar and the radio- 
sonde computations,. Figure 10 shows 
the two heights as a function of time and 
the differences plotted to the right. In 
this example the height accuracy of the 
radiosonde is quite good from the sur- 
face to 25 km. Above this height the 
errors increase rapidly to about 
1 km at 35 km. Errors 2 to 3 times 
this large were found in some ascents 
and these were attributable to errors 
in pressure. A more complete analysis 
of the height errors is possible, time permitting. 


HEIGHT COMPARISON 24/7/79 3RD ASCENT 



INDTE:. RflDflft H7 OFFSET 0T J00HJ 


Figure 10. Error in geopotential heights 
computed from radiosonde. 


An example of the filtered profiles is presented in Figure 11, left and right pages. 
Starting with the upper left side, the profiles from the low-pass Filter No. 3 are plotted 
across the page. The profiles include Relative Humidity, Temperature, Zonal Wind, 
u component and Meridional Wind, v component. Winds from the west and south are 
positive in these plots. Along the second row on the left, profiles of the 0. 9-2. 7 km 
band-pass are plotted against the same height scale, but with appropriate changes in 
the scales on the abscissa. At the lower left, the band pass components are added to 
the Filter No. 3 profiles to show the relative amplitudes and how the waves add 
structure to the profiles. 
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0.5 MIN FIlTfW APPLIED TO MANAMA DATA 


7/2*/i9/7 


Z(M> 

UT 

l (SFC > 

0ELZ 

pp 

tkmp 

KH 

h A V 

TV 

THV 

TnVL 

MD 

MS 

u 

V 

. M 
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H*S 

oPOT/M 

GPOT/M 
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(C) 


(Mm) (6/KG) 

(ft) 

(ft) 

(K) 
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(M/S) 

( M/5 ) 

(M/b) 

(M/b) 


0 

16 

?9 

0 

0.00 

0.00 
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28.20 

83.00 

32.511 

20.366 

305.10 

304.6b 
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3*0.00 

2.00 

.68 

-1.88 


1 

16 

29 

35 

99.82 

-.00 

993.96 

25.70 

92.37 

30.313 

19.319 

302.39 
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354.26 






2 

16 

29 

58 

199.61 

.01 

982.82 

25.18 
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80,851 

19.565 
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303. 4U 

365.48 
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6,04 

.07 
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4.31 

3 

16 

30 

21 

299.41 

.01 

971.77 

24.2/ 

97.99 

29. 723 

19.249 

300.93 

303.40 
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857.7b 

5.73 

.22 

-5. 73 

4,40 

* 

16 

30 

43 

399.22 

.01 

960.82 

24.11 
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27.003 

1 7.668 

300.48 
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361.27 

859. bb 

5.95 

.04 
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4,60 

b 

lb 

31 

5 
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.00 

949.97 

23.54 

68. 78 

2b. 769 

17.047 

299. 78 

304.22 

349.93 

859.45 

5.92 

.06 

-5.91 

4,3b 

6 

16 

31 

29 

598.80 

.02 

939.23 

22.94 

88.04 

£4 . 6*6 

16. *85 

299.08 

304.49 

348. 7* 

880. bi 

6. 74 

.94 

-5.66 

4,22 

7 

16 

31 

52 

69b. b0 

.01 

928.57 

22.68 

77.04 

21.222 

14.3*0 

298.42 

304.82 

343.67 

3*1.2* 

6.75 

1.65 

-5,44 

4.29 

H 

16 

32 

16 

/9H.38 

.01 

918.03 

22.42 

73.36 

19.903 

13.397 

298.03 

305.41 

342.44 

3*4.02 

6.08 

1.67 

-5.84 

4.J2 

9 

16 

32 

38 

898. 19 

-.00 

907.58 

21.57 

80 . 95 
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14. *16 

297.32 

305.69 

3*4.94 

8*5. *b 

b.01 

1.51 

-5.82 

4.32 

10 

16 

33 

2 

997.92 

.05 

897.23 

20.84 

79.88 

19.668 

1 3.749 

296.46 

305.80 

343.31 

3*2.1* 

5,78 

1.77 

-5.50 

4.22 

11 

16 

33 

26 

1097.38 

.36 

687.00 

20.02 

65.52 

20.020 

14.160 

298. 72 

30b. 03 

344.63 

3*0, *5 

5.71 

1.34 

-5.55 

4.13 

12 

16 

33 

b» 

1195.19 

2.33 

877.03 

19.41 

Hb.95 

19.399 

14.018 

295.07 

306.35 

3*4.0 7 

850. *b 

5.08 

.31 

•5.02 

4.06 

13 

16 

34 

15 

1298.70 

-1.41 

666.58 

19.22 

79.68 

17. 74* 

1 2 .635 

294.67 

306.98 

342.26 

7.2b 

4.82 

-.61 

-4.76 

4.05 

14 

16 

34 

40 

1401 .64 

-4.77 

856.26 

19.17 

7d.*5 

17.424 

12.735 

294.61 

307.98 

343.20 

21.27 

4»b9 

-1.70 

-4.37 

3.67 

IS 

16 

35 

10 

1499.96 

-3.13 

846.59 

18.54 

73. 30 

15.651 

11.580 

293. 76 

308.09 

340 • 22 

87. b2 

4.12 

-2.52 

-3.27 

3.47 

16 

16 

35 

37 

16U2.39 

-5.79 

836.56 

18.22 

72. bl 

lb. 19 7 
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*b . 06 

3,95 

-2.94 
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3.62 

17 

16 

36 

5 

1 706.46 

-10.10 

626.48 

17.76 

72.76 

1*. 796 
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309.32 

3*0.64 

61.19 

*.*7 

-3.92 

-2.15 

8.97 

IB 

16 

3b 

28 

1 798.95 

-2.83 

617.61 

17.12 

7* .68 

14.573 
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6. 01 

-6.b* 

-1.44 

4.45 

19 

16 

36 

50 
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-2.65 
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15.126 

11. 725 

291.74 
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7.97 

-7.95 

-.59 

4*Od 

20 

16 

37 

11 
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16.24 
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14.092 
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9.78 

-9.76 

• 48 

4.77 

21 

16 

37 

32 
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-7.71 

769.00 
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12. /93 

10.148 
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-10.36 

.27 

4.77 

22 

16 

3/ 

52 
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87,9b 

10. 80 
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-.36 
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23 

16 

38 

14 
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-2.10 

771.17 
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76.34 

12.130 

9.842 

288.80 

311.08 

338.96 

89.51 

10.97 

-10.97 

-.09 

4.75 

24 

16 

38 

35 

2406.76 

-12.17 

761.21 

13.65 

64.30 

10.0*9 

8,233 

288.26 

311.65 

335.24 

92.9b 

10.72 

-10.70 

.56 

4.7b 

2S 

16 

38 

55 

2509.66 

-15.35 
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13.54 
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7.938 
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312.56 

335.44 
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11. bb 

-11.85 
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4.73 

26 

16 

39 
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-4.03 
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39.39 

6.96* 
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89,95 
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4.69 

27 

16 

39 

38 
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7 35. 14 
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4.54 

28 

16 

40 

1 
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34.09 
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-13.64 
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4.40 

29 

16 

40 

24 
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-11.67 
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31 3.8/ 

33 3.15 
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-13.64 

-1.50 

4.45 

30 

16 

40 

46 
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-12.92 
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10.37 
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6. d23 

284. 71 

314.16 

334.06 

82.71 

13.31 

-13.20 

-1.69 

4.49 

31 

16 

41 

8 
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-9.69 
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9.79 
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6.688 

284.14 

314.58 

334.69 

b2 . * 1 

12.00 

-11.69 

-1.58 

4.53 

32 

16 

41 

30 
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-10.84 

692.17 

9.08 
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6.544 

283.3/ 
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334.00 
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11.3/ 

-11.24 

-1.75 

4.59 

33 

lb 

41 

52 
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-11.65 
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8.26 
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7.629 

6.968 
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335. *3 

7b. 89 

11.41 

-11.18 

-2.30 

4.56 

34 

16 

42 

14 

3403.68 

-11.88 
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7.69 
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6 ,306 
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281.88 

315.33 
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80.19 
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-11.56 

-2.00 

4.53 

3b 

16 

42 

36 

3506.93 
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7.07 
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281.13 

315.05 
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12.30 

-12.12 

-2.12 

4.59 

36 

16 

42 

57 

3610.92 

-19.71 
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3.91 1 

2d0 . 33 
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327.97 

82.56 

11.2b 

-11.19 

-1.46 

4.73 

37 

16 

43 

18 

3708.34 

-17.44 
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6.74 
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3.521 
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9.27 
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4.72 

38 

16 

43 

40 
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-15.56 
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35,21 
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.55 

4.62 

39 

16 

44 

1 
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-15.95 
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7 • bb 
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1.97 
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40 

16 

44 

23 
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4.60 
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4.b8 

41 

16 

44 

44 
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4.858 
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318.56 
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1 lb. *9 
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-6.92 
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4.70 

42 

16 

4b 

5 
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-14.55 
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b.18 
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4,50 

4.83 

43 

16 

45 

26 
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-18.82 
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2.81 
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5.00 

44 

16 

4b 

45 
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3.79 
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45 

16 

46 

5 
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3.893 
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332.63 
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7.87 
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5.10 

*6 

16 

46 

25 
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320. b7 
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7.92 

- 7.3u 
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16 

46 

44 
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274.00 
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7.98 
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3.27 
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4b 

16 

47 
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568.49 
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16 

47 
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3 MIN FIl'TF* APPLIED TO PANAMA DATA 7/24/1977 
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WH 

t 
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WO 
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0 

16 

?9 
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0.00 
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29 
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2 

16 

29 

55 
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3.82 

3 

16 

30 

20 
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-5.33 

4.05 

4 

16 

30 
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6.25 

-.57 
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4.28 

5 

16 

31 

7 
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The sums, the light lines in the lower left row, are repeated in the upper right row. 

These profiles represent the low pass for the 0. 9 km filter. Filter No. 2„ Below them 
are the 0. 15-0. 9 km band-pass profiles and in the lowest row are the profiles resulting 
from applying Filter No. 1 to the unfiltered data. Note that the 0„ 3 km heading for the 
middle row and the 0. 3 km (one minute) heading for the lower row are incorrect by a 
factor of two. Filter No. 1 has a 50% response at 0. 15 km which is equivalent to 
0. 5 minute for an average balloon ascent rate. 

These are the basic diagrams for displaying the results. They are completely automated 
and were drawn to scale for ease in reproduction. All ordinates are identical and 
correspond to geometric not geopotential heights. Each figure contains a vast amount 
of information and the compressed vertical scale obscures some of the phase relations, 
especially those of the shorter wavelength components. But by being able to see the 
entire set, without turning pages and losing one's visual reference, it has definite 
advantages for those who "see" what they are looking at. 

The band-pass profiles clearly establish the presence of waves with a short vertical 
wavelength. A survey of all of these plots. Figures 15 to 54, indicates that they are 
always present and that they contribute to the fine structure of the temperature and 
relative humidity profiles as well as to the u and v velocity components. In general, 
amplitudes of the two velocity components are comparable within a given band pass, 
with decreasing amplitudes associated with decreasing wavelengths. There is evidence 
also for nodes at approximately 11 and 22 km in these band-pass profiles, the latter 
being the approximate height of the easterly jet. 

Evidence of left and right circular polarizations, elliptical and linear polarizations are 
to be seen in the 0=,9-2„ 7 band-pass profiles. For example, in Figure 11, left side, 
middle row, the band -pass wind at approximately 1 km is from the west at 2. 5m sec - ^. 
With increasing height the wind component changes to blow from the south, then from 
the east, then from the north and finally back again from the west at approximately 4 km. 
This counterclockwise rotation (cyclonic), by analogy to the righthand rule, will be 
called right circular polarization,, At approximately 15 km the polarization is reversed 
and at approximately 28 km the u and v components are almost linearly polarized. 

Internal gravity waves are linearly polarized, but only where the larger scale winds 
have little vertical shear, i„e. , where the thermal field is quasibarotopic. As shown 
by Danielsen (1966), vertical advection of the larger scale wind shear produces a 
relative phase shift between u and v which approaches quadrature when the shear is 
large. 

Amplitudes of the 0. 9 to 2. 7 km band-pass components are often 3 to 5 m sec~l, definitely 
nontrivial. Associated with these wind oscillations are temperature perturbations of 
approximately 0. 5° in the troposphere which increase to approximately 1 to 2° in the 
stratosphere. Apparently the adiabatic compression and expansion, not horizontal 
advection, dominates for the local temperature change. Note that the relative humidity 
changes are quite the reverse. They maximize in the troposphere and indicate a 
dependence on horizontal gradients. 
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The same pattern of changes in perturbation amplitudes versus height can be found in the 
shorter wavelength band-pass, middle row, right figure. Again the tendency for a node 
at approximately 22 km is evident. However, all the amplitudes are about a factor of 
2 smaller except those of relative humidity. This difference supports the relative im- 
portance of horizontal gradients for humidity perturbations. It is highly probable that 
these waves and those of the 0. 9 to 2. 7 km band-pass have high frequencies and con- 
sequently cannot be resolved from synoptic radiosonde, rawinsonde observations. 
However, it is clear that they contribute a complexity to tropopause identification. The 
low-pass profiles in the upper left are more representative of the "synoptic" scale 
structures, but the temperature profile is curved rather than discontinuous in slope, 
i. e. , the tropopause is a second order discontinuity in T or a place where the curvature 
increases rapidly. 

It is evident from the band-pass profiles that the amplitudes of the short, vertical wave- 
length oscillations are rather large. Velocity amplitudes of 2 to 5 m sec - ^ have a 
potential for affecting the transport of trace gases and aerosols between the troposphere 
and stratosphere and within the stratosphere itself. The relative importance of these 
waves to transport at tropical latitudes remains to be established, but their presence 
can no longer be neglected nor denied. A significant eddy-type transport in the tropical 
stratosphere was deduced by Louis (1974) as being essential to maintaining the ozone 
concentrations below the chemically active altitudes against the upward and difluent 
transports by the mean circulations. 

Having demonstrated the presence of the short, vertical waves, we decided to apply two 
additional low-pass filters to separate the longer vertical, wavelength components. 

Despite the broad transition zone of these two filters, especially that of Filter 5, 

Figure 12 reveals some interesting results. The 24-8 km band-pass, middle row, left, 
shows a distinctly different wavelength for the zonal than for the meridional components. 
Oscillating between m sec'^, the zonal winds have a vertical wavelength of 
approximately 16 km. A similar wavelength dominates the temperature deviations, 
with the phase lag similar to that of a Kelvin wave (Wallace and Kousky, 1968). The 
linearly-polarized Kelvin wave with velocity components restricted to the X, Z plane 
is a special type of internal gravity wave restricted to equatorial latitudes where 
f = 2uu sin<t>, the vertical component of the earth's vorticity, approaches zero, w being 
the angular speed of the earth's rotation. In the Kelvin wave, these zonal winds are in 
geostrop hie balance. Although f approaches 0 as $ approaches 0 , Ug remains finite 
when the pressure surfaces have a maximum ( or minimum) elevation at the equator. The 
transition from westerly to easterly winds then requires negative temperature departures, 
as seen in Figure 12, which changes the curvature of the pressure surfaces from negative 
to positive with increasing heights. 

Because the 16 km wave is linearly polarized in the zonal direction with temperature 
departures in quadrature, one is tempted to identify it as a large vertical -wavelength 
Kelvin wave. However, any such temptation is thwarted when one examines its local 
phase propagation. To qualify as a Kelvin wave the local vertical phase propagation 
should be negative or downward and the horizontal phase propagation should be positive 
or eastward. But this 16 km wave, although fluctuating diurnally up and down, has only 
a very small downward phase propagation. For example, Figure 13 shows the same 
wave eight days earlier in almost the same phase as in Figure 12. We will establish 
its probable identity later. 
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In contrast to the quasistationary phase of the zonally polarized 16 km wave, the meridional 
components of an approximately 11 km wave have exactly reversed phase during the eight 
day period. A similar analysis for the 19th of July, not shown, indicates that this 11 km 
wave is produced by interference between upward and downward propagating wave modes. 

In other words, wave reflection or ducting is apparently taking place. The amplitude of 
this north-south oscillation is smaller by a factor of 3 or 4 than that of the east-west 
oscillation, but again, it is significant. 

In addition to these two large-scale waves, there is evidence in the middle row, right 
of Figures 12 and 13, for rather large amplitude velocity components associated with 
approximately 5 kra wavelengths. This vertical wavelength approximates that required 
for the Kelvin and gravity -modified Rossby waves but the phase of the temperature 
deviations doesn't fit well to either wave mode. In any case, the amplitudes of the 
wind components (less than or equal to 5 m sec - -*-) are quite large, therefore, their 
potentials for transport cannot be ignored. 

From these data it is an easy matter to compute the temperature or potential temperature 
flux in the zonal and meridional directions. The profiles in Figures 11, 12 and 13 
suggest that these fluxes will be oscillatory but certainly significantly large in magni- 
tude, with the largest fluxes produced by the largest vertical waves. It is interesting 
to note that the amplitude of the velocity deviations in meters per second is approximately 
equal to the vertical wavelength in kilometers; also, as is typical of internal 
gravity waves, the temperature perturbation in degrees is about 50% of the velocity 
perturbation in meters per second. From these general estimates one can deduce that 
the temperature flux increases as L^, where L is the vertical wavelength. Since 
we have evidence here of waves with L ranging from 1 to 16 km, the heat flux pro- 
duced by the largest waves are 100 to 2b0 times those of the 1 km waves. 

Estimates of the meridional diffusion coefficient, Kyy a diagonal component of the 
diffusion tensor, can be made also from these data. For example, if we scale K by r 


K — - u'v't (11) 

yy v ’ 

use v' = 10 m sec -1 , and following Danielsen and Deaven (1974), we set t = 4xl0 4 sec we 
obtain a Kyy— -10^ cm.^ sec - *. This agrees with the magnitude deduced by Louis (1974) 
for the equatorial zone. 

Returning now to the enigmatic 16 km zonal wave which is quasistationary locally, it 
was recognized that, except for the forcing function, the generalized wave equation for 
stationary mountain waves in a nonuniform stability and wind regime (Danielsen and 
Bleck, 1970) would be applicable. This equation determines the total wave number 
squared as a complicated function of the stability and its vertical derivative, and the 
wind speed and its vertical derivatives. Then, knowing x (Z) and k£ one can solve 
for kg (Z). The computed k values, based on the 24 km filters, yields k|r values 
which are about 10 times too large, i. e. , the L^'s would be closer to 5 than to 16 km. 
Again, the more general solution does not permit a free mode with such a large 
vertical wavelength. 
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However, this wave, which dominates the zonal wind profiles obtained from radar track- 
ing of the radiosonde balloon, also dominates the wind profiles determined from radar 
tracking of the rockets which were released and tracked from the same site. Figure 
14 shows the zonal wind profiles from our 8 and 24 km filters for the 24th of July 
with the rocket-derived winds from Schmidlin and Kloos (1978). Note that the latter 
and the former agree with respect to phase, but the rocket winds have a smaller 
amplitude. The mean profile, in the center, was drawn by hand and then these values 
were subtracted to produce the deviation profiles, at the right. 


ZONAL WIND PROFILES 24JULY 1977 



Figure i4. Zonal wind profiles from 24-8 km band-pass and from radar tracking 

of rocket. 


From Figure .14, right, it is apparent that this oscillation extends from the surface through 
the troposphere, stratosphere, and into the mesosphere. The vertical wavelength varies 
from approximately 14 to 18 km and the amplitude is smallest near 40 km. Clearly, 
it dominates the zonal wind profile. After examining the entire set of rocket winds 
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it was obvious that the miniraae and maximae fluctuated slightly in altitude and speed 
but there was no significant downward propagation. Thus, the rocket and radiosonde 
winds agree in all respects for this approximately equal -to-16 km wave. 

The obvious implications are that this wave is forced and has a very long period. In 
general, it seems to fit the properties of the semiannual waves analyzed by Wallace 
(1973). Indeed, the profile for the summer of 1963 in Wallace's article is extremely 
similar to that in Figure 14. The mesospheric and upper stratospheric easterlies 
agree also with those illustrated in Figure 9. 6 of Danielsen and Louis (1977). Here, 
the easterlies are attributed to the divergence produced by the heated hemisphere. 

The easterlies propagate across the equator, then reverse to become westerlies in the 
cooling winter hemisphere. Therefore, tropical easterlies are produced by both the 
northern and southern hemispheric summers. While they are approaching the equator 
in the heated hemisphere the winds from the east are approximately balanced by a 
warm anticyclone in the pressure field. However, a small ageostrophic meridional 
component advects the easterly winds across the equator where they experience a 
reversal of the coriolis force. As this force increases with increasing latitude, the 
resulting imbalance of forces produces rapid anticyclonic turning which reverses the 
flow's direction to westerly. Then, this westerly flow is approximately balanced by 
the cold cyclone in the pressure field of the winter hemisphere and the small ageostrophic 
meridional component continues the convergence toward the cooling pole, with an 
associated increase in zonal wind speeds. However, although this easterly momen- 
tum source during the solstices may explain the semiannual oscillation above about 
35 km, it does not explain how the vertical wavelength is established. 

SUMMARY 

High resolution wind, temperature and relative humidity profiles have been derived by 
objective computer methods from a unique set of observations made in the vicinity of 
the Intertropical Convergence Zone. The set is unique because: 

1) radiosonde-ozonesonde observations were made approximately every six 
hours for the latter half of July, 

2) the balloon-borne system was simultaneously tracked by the conventional 
GMD-1 and by a NIKE radar, and 

3) all data were digitally recorded at a rate of one or two per second. 

Taking full advantage of the rapid recording rate, the data were objectively filtered by 
applying a series of non- recursive low-pass filters. The filtered data are presented 
as vertical profiles of the east-west and north-south wind components, the temperature 
and relative humidity. 

The first filter removes the very short vertical waves and the majority of the instru- 
mental and recording errors. The remaining filters remove progressively longer and 
longer waves to reveal the large-scale structures of the tropical atmosphere. The 
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largest scale wave, with a vertical wavelength of 16 *2 km, dominates the zonal winds 
and the temperature profiles from the earth's surface through the troposphere, strato- 
sphere and well into the mesosphere. This wave, with some of the properties of a Kelvin 
wave, does not have the proper intrinsic or local phase propagation for a Kelvin wave. 

It probably represents a semiannual oscillation generated in the upper atmosphere by 
heating in the summer hemisphere. 

Of special significance to the initial research objectives are the profiles of the band -pass 
components, which firmly establish the existence of a broad range of mesoscale internal 
waves. The amplitudes of these waves decrease with decreasing wavelengths, but at 
all wavelengths they are significant. In particular, mesoscale waves in the wind with 
amplitudes greater than 5 m sec - '*’ are evident, especially in the stratosphere. The 
importance of these waves to stratospheric-tropospheric exchange and to eddy trans- 
ports remains to be determined, but their potentials are now established. 
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